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The purpose of this study is to determine the experimental uncertainties in measuring mass and cross section param- 
eters of SUSY particles at a 500 GeV Linear Collider. In this study SUSY is a point in the focus point region of 
mSUGRA parameter space that is compatible with WMAP constraints on dark matter relic density. At this study 
point the masses of the squarks and sleptons are very heavy, and the only SUSY particles accessible at the Linear 
Collider would be the three lightest neutralinos, and the two lightest charginos: Xi> X21 X31 X\ 1 X2 > wnere X? i g the 
lightest supersymmetric particle (LSP). The charginos or neutralinos may be pair produced, and the subsequent decay 
cascades to the LSP allow us to measure the SUSY couplings and mass spectrum. We find that by looking for the 
qq i signature 2 jets plus 2 leptons plus missing energy we can determine the mass of the LSP to within 1 GeV uncertainty 

and that the mass differences of \% an d X3 with the LSP mass can be determined to better than 0.5 GeV. 

o 
o 

O 1. INTRODUCTION 

<D 
<D . 

This study represents work in an effort to establish a connection between measurements at a future linear collider 
and cosmological data. We assume that the LHC has run and discovered SUSY and that a 500 GeV linear collider 
has run to make precision measurements. The ultimate goal of this study is to determine whether or not the Linear 
Collider would be able to make precise enough measurements of mass and cross sections of SUSY particles in order 
to determine if the relic density of the LSP would be consistent with WMAP constraints [l| , and that the SUSY 
LSP actually is the dark matter. 

We have chosen the scenario where the discovered SUSY is mSUGRA in a focus point region with Mi /a = 300 
GeV, A — 0, tan/3 = 10, \i > 0, m t = 175 GeV, m — 3280 GeV. In this focus point region the masses of the 
squarks and sleptons are very heavy, and the only SUSY particles accessible at the Linear Collider would be the 
three lightest neutralinos, and the two lightest charginos: Xit X,2> Xssj Xii xt > where neutralino Xi is the LSP. The 
signed mass parameters and the production cross sections of the charginos and neutralinos at this focus point 
region are given in tables U and[n] Note that the mass parameters are signed. The sign of the mass parameter will 
affect the decay distributions but not the kinematics. 

This paper represents preliminary results on the expected experimental uncertainties in determining the mass 
parameters and production cross sections of SUSY particles at the linear collider. Specifically, we will discuss the 2 
jet plus 2 lepton plus missing energy signature of X2X3 production. 
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"This talk was given by R. Gray, describing ongoing work performed in collaboration with the other authors. 
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Table I: Neutralino and Chargino Mass Parameters 
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Table II: Production Cross Sections of SUSY particle pairs in collisions of 95% left or right polarized electrons on unpolarized 
positrons. 
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2. SIMULATION AND ANALYSIS CUTS 
2.1. SUSY Signature 

We will analyze the production of X2X3- Since the squarks and sleptons are so heavy at this point in the parameter 
space, the dominant means for the \2 an d X3 to decay is a transition to Xi (the LSP) by radiating a Z. Since the 
mass differences between x", X§ an d X? are all smaller than the Z mass, the Z is virtual. In a detector, the two 
LSPs produced in the decay of the X2X3 pair w iU be seen only as missing energy. The two Zs produced in the 
decay will each create one of the following: 11, vv, qq. A 11 pair will be visible as tracks, a vv will contribute to the 
missing energy, and a qq will create at least 2 jets. Therefore, we will be looking for one of the following detector 
signatures for X2X3 production: 2 Jet + -E m i ss , 2Lepton + -E^igg, 4 Jet + -E m j ss , 4Lepton + -E m i ss , or 
2Lepton + 2 Jet + £^33 • For X2X3 production we will use the mode 2Lepton + 2 Jet + -E^igg because it has the 
smallest backgrounds from the other SUSY and Standard Model processes. 



2.2. Event Generation 

SUSY events were generated at Cornell using Isajet 7.69 This included production via e + e~ collisions of the 
following neutralino and chargino pairs: xtxi 1 Xi~X2~j X1X3; X°X3- The Standard Model backgrounds for this 



X2 i X1X3 
1.22 Q. 



analysis were generated at SLAC using WHiZaRD 1.22 4]. The background include contributions from the 
following collisions: e + e~, e~7, e + 7, 77. The dominant backgrounds for this analysis generated by these collisions 
are tt, W + W~ and ZZ. The beam in these event generators is a ±95% polarized electron beam on an unpolarized 
positron beam. The generators include effects from initial state radiation (Beamstrahlung and Bremsstrahlung) 
which can reduce the energy available for SUSY production. 

2.3. Detector Simulation and Reconstruction 

The detector simulation was done using LCD Root FastMC, SD MarOl, which is a root based fast Monte Carlo 
program produced at SLAC. Once we have tracks and showers the first task is to identify leptons that are isolated 
from jets and are candidates for having come directly from a neutralino or chargino decay. It is assumed that at the 
linear collider identification of electrons and muons will be near perfect and so lepton ID uses the Monte Carlo 
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truth information. These leptons are candidates for direct production from a SUSY decay if the lepton momentum 
is greater than 10 GeV and the total energy from other particles within 20 degrees of the lepton is less than 2 GeV. 
After these isolated candidate leptons are identified, perfect track-shower matching is employed to approximate the 
"energy flow" algorithms that would be used at a real detector and the remaining 4-vectors are divided into jets 
using the Durham algorithm with a y cu t value of 0.004. 

We will use b-tagging in this analysis to help identify tt backgrounds. Since the detector simulation used in this 
study is not advanced enough for an honest b-tag, Monte Carlo truth information was used. The jets found by the 
Durham algorithm were matched to the generator level jets by an energy weighted comparison of the lists of stable 
particles in the "found jet" and the list of particles in the generator level jet. A match quality is defined by 
equation ^ where Ei is the energy of particles in list i, Ei a j is the energy of particles that are listed in both list i 
and list j, and Ei n j is the energy of particles listed in list i but not in list j. If the match quality is 1, the two lists 
are the same, if the match quality is —1 they have nothing in common. The generator level jet with the highest 
match quality for a "found jet" is tagged to that found jet. If the "found jet" is matched to a jet that came from a 
b quark, then we say it has a 50% chance of being b-tagged. 

IE ■ — E- ■ — E 

Qij ~ E i +E j W 

2.4. Analysis Cuts for xlxl Ifaqxlxl 

The analysis cuts for the 2 jet plus 2 lepton plus missing energy signal are as follows. There must be two leptons of 
opposite charge and same flavor that both passed our isolation cuts. There must be 2 or 3 jets found in addition to 
the isolated leptons. The missing energy must be greater than 275GeV. The transverse momentum must be greater 
than lOGeV. The cosine of each track and jet with the beam axis must be less than 0.95. The energy of each jet 
and letpon must be less than HOGeV. Each jet is then anti-btagged, or the event is cut if either jet appears to have 
come from a b quark. This anti-btag cut removes backgrounds from tt. 



3. MASS PARAMETER AND CROSS SECTION MEASUREMENT 

3.1. Cross Section Measurement 

The cross section will simply be measured by counting the number of events that pass our cuts with backgrounds 
subtracted, divided by the efficiency for passing those cuts, and divided by the luminosity. With 500 fb _1 of 
luminosity divided evenly between left and right polarized beams we find statistical uncertainties on polarized cross 
sections of 4% and 5% respectively. 

3.2. Mass Parameters from Invariant Mass Distribution 

The mass difference between an d X?> an d the mass difference between x§ and xi, may be determined by fitting 
the invariant mass distribution of the two leptons in the 2 letpon plus 2 or 3 jet data. For any decay Xi — ► UXf the 
invariant mass of the two leptons will have a maximum possible value of the mass difference between Xi and Xf. 
Since x§ an d X§ have different masses we will see two edges in the mass distribution. One edge will be for the \\ 
and the other for the \%- I n order to use the whole distribution, rather than just the edges, in our mass 
determination, we will fit all the data to the expected invariant mass distribution. This expected mass distribution 
will have very weak dependence on the LSP (x?) mass, but will depend on the relative sign of the \\ an( i X3 mass 
parameters with the x\ mass parameter. You can see from table [i] that x\ and x% mass parameters have different 
relative signs. We find that the efficiency for measuring the di-lepton is mostly flat in invariant mass beyond 10 
GeV, so we will make no efficiency corrections to the distribution prior to performing the fit. 
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2J2L, 3 Sign Possibilities Shown 
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Figure 1: This shows the invariant mass distribution of combined 250fb~ 1 left polarized and 250 fb _1 right polarized data. 
The magenta is the signal. The yellow is standard model backgrounds, and the blue is background from chargino production. 
The three fits are for the different sign assumptions for the neutralino mass parameters. The red line is for (-,-), the blue line 
is for (+,+), and the black line is for (+,-)• The correct sign assumption is (+,-). 



The resulting invariant mass distribution can be seen in figure ^ We expect two edges in the distribution. We 
expect one at 58.6 GeV for X21 an d one a * 82.3 GeV for the X3. Both edges are visible. 

The mass distribution of the background is fit to an order 2 polynomial. The fit to the background distribution is 

combined with the theoretical prediction for the SUSY component and an unbinncd log likelihood fit is performed. 

The fit is done three times for the three possible combinations of relative signs of X2X3 to the Xi- (t)j (+>+)> 

(+,-). We find that the possible sign combinations are easily distinguished in the data, and that the best fit for the 

correct sign combination, (+,-), is at least 13 standard deviations away from the best fit of the incorrect sign 

combinations. The fits for the different sign combinations are shown in figure ^ 

The log likelihood fit for the mass differences yields M^o — M^o = 58.8GeV ± 0.3 and 

M^o — M^o = %2.2hGeV ± 0.2. The generator level values for these masses are 58.6 GeV, and 82.3 GeV. 

ft should be noted that theoretically we could do the same study with the invariant mass of the hadronic Z decay. 

However, reconstruction of di-lepton mass distribution will have much less uncertainty than reconstruction of the 

di-jet mass distribution. Also the efficiency and sculpting of the di-jet distribution is much more complicated than 

for the di-lepton due to the jet algorithms. 
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3.3. LSP Mass From Energy Distribution 



The invariant mass distribution of the two leptons had very weak dependence on the LSP mass. However, the 
distribution of the energy of the two leptons will depend on the LSP mass. It will also depend on the X2 an d X3 
masses, the theoretical invariant mass distribution of the two leptons, the beam energy, the spectrum of initial state 
radiation, and the spin correlations in the neutralino decays. 

The energy distribution is numerically calculated by integrating over all of the above mentioned distributions. It 
should be noted that the version of Isajet used in this analysis docs not calculate the spin correlations of the 
neutralino decays Q ■ The spin correlations of the neutralino decays affect the angular distribution of the virtual Z 
in the rest frame of the parent neutralino along the boost direction. This angular distribution is flat in the Monte 
Carlo because of the lack of neutralino spin correlations. This translates to a flat energy distribution in the lab 
frame for a given value of di-lepton invaraint mass and beam energy. 

The data and the calculated energy distribution are shown in figure [21 The standard model background 
distribution is fit to an order 3 polynomial and is combined with the calculated signal distribution for a log 
likelihood fit to the data. In the fit the x? (LSP) mass is allowed to vary, and the mass differences of the xt an d the 
%3 with the LSP mass are allowed to vary within 3 standard deviations of the measured mass differences obtained 
from the invariant mass distribution fit. The fit finds an LSP mass of 108.3 GeV ± 1.0 GeV. The generator level 
value of the LSP mass is 107.7 GeV. 



4. Conclusions and Future Directions 

We have presented preliminary results for the expected uncertainties that a 500 GeV linear collider would have in 
measuring mass and cross sections from X2X3 production in a focus point region using the 2 lepton plus 2 or 3 jet 
plus missing energy detector signature. We have found that for 250 fb _1 of left polarized data and 250 fb~ of right 
polarized data we could determine the cross sections to within a uncertainty of 4% and 5% respectively. By using 
the invariant mass distribution of the two leptons of the combined 500 fb -1 data set we could measure the mass 
differences M^o — M^o and M^o — M^o with 0.3 GeV and 0.2 GeV uncertainty respectively. By using the energy 
distribution of the two leptons we can measure the Xi (LSP) mass with an uncertainty of 1.0 GeV. 
Other work currently being done by the Cornell/Florida LC Cosmology Group is the study of X1X2 production, 
X1X1 production, and X1X2 ■ There is also work being done to convert the uncertainties in mass and cross section 
measurements into uncertainties in the relic density. If the uncertainty in the relic density of the LSP is small 
enough, we can compare our prediction to WMAP constraints and determine whether or not the LSP actually is 
the Dark Matter. 
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2J2L, Di-Lepton Energy, After cuts, 500fb 
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Figure 2: This shows the energy distribution of the two leptons in the 2 lepton 2 or 3 Jet and missing energy decay signature. 
The magenta color is the signal, the yellow is the contribution from standard model backgrounds and the blue is background 
from chargino production. The curve on the plot is the calculated energy distribution for the signal. 
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